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Abstract 

We characterize group symmetries of poly-phase complementary code matrices (CCMs), which we use 
to classify CCMs in terms of their equivalence classes. We also present classification results for CCMs of 
dimension N x 4 where N = 2,3, 4, 5,6. Finally, we present a new construction to generate quad-phase 
CCMs from ternary CCMs and compare this to other existing constructions that focus on generating 
CCMs from those of smaller dimensions. 

1 Introduction 

In phase-coded pulse-compression radar, signal waveforms can be represented by a string of values (called a 
code) consisting of roots of unity. In order to correlate incoming with outgoing signals these codes should 
have nice sharp peak-sidelobe characteristics in terms of their auto-correlation functions. Usually in radar 
and communications, Barker codes and Golay pairs of codes are employed, but there is increased interest 
in more general codes called complementary code sets, i.e., complementary sets of codes whose composite 
auto-correlation function has zero sidelobe levels mm- When expressed in matrix form, complementary 
code sets are referred to as complementary code matrices (CCMs). 

In this paper we characterize known symmetries of CCMs in terms of their relations. Previously, Golomb 
and Win investigated symmetries of a single polyphase sequence. These symmetries were later extended 
to CCMs by Coxson and Haloupek [T]. Let M be a CCM. Here are the five known symmetries of M that 
preserves its CCM-property: 

(i) Column multiplication by a unimodular complex number. 

(ii) Column conjugate reversal. 

(iii) Matrix conjugation. 

(iv) Progressive multiplication by consecutive powers of a unimodular complex number. 

(v) Column permutation. 

When viewed as group generators these five symmetries are non-commutative in general, e.g., column mul¬ 
tiplication does not commute with matrix conjugation. Therefore, it is important to characterize their 
relations, which we use to classify CCM’s in terms of their equivalence classes. 
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Our results describing the group relations between the five symmetries above extend those of Coxson [7] 
who determined the group structure for symmetries of Barker codes. Moreover, we obtain an upper bound 
on the corresponding group generated by these symmetries, which we call the complementary group G. For 
p-phase N X K CCMs, we establish in Section [3| (Theorem |3.3[ ) that 

|G| < 2^+^p^+^Kl 

In Section we extend Coxson and Russo’s |2] efficient exhaustive search algorithm for binary CCMs 
to p-phase CCMs. This algorithm was implemented for quad-phase x 4 CCMs in Section to obtain 
a classification of all equivalence classes for N = 2,3,4, 5, 6 (see Table for a description of the number 
of equivalence classes and those that are represented by Hadamard matrices). This extends Gibson’s |3] 
classification results for quaternary Golay sequence pairs. In the same section, we present a method to 
construct quad-phase CCMs from dual-pairs of ternary CCMs (dual in the sense that their commutator 
is diagonally regular). Lastly, we present in Table a list of the number of equivalence classes whose 
representatives are dual-pairs. 


2 Preliminaries 

In this section we present definitions of complementary code matrices based on the auto and row-correlation 
functions. Let M he a N x K p-phase matrix, i.e., one whose entries consists of p-th roots of unity. At times 
we will represent M in several different ways: coordinate-wise with M = (l<n<A^, 1 < k < K), 

column-wise with M = [xi,X2, X3, ..., xk], or row-wise with M = [ri, r 2 , ^ 3 ,... ■ 

Definition 2.1. The aperiodic autocorrelation function of a code x = (ui,... ,ajv) is defined to be 

{ N-j 

aia^+j, if 0 < j < TV - 1 ; 

( 1 ) 

A.(-j), if - iV + 1 < j < 0. 

Note that Aa;(0) = \x\‘^ = N. Also, a desirable autocorrelation function should have a sharp center peak 
(j = 0 ) and low sidelobes (j 7 ^ 0 ). 

Example 2.1. Let x = { —1,—1,—i}. Then 

Aa;(-2) = -l(z) = -i 

A,(-l) = -l(-l)-l(z) = l-f 
A,(0) = -l(-l) - 1(-1) - i{-i) = 3 
A,(l) = -l(-l)-l(-f) = l+i 
A,(2) = -l(-f) = L 

Example 2.2. An example of a Barker Sequence, where the | Ax{j) |< 1 for all j 7 ^ 0, is the following: 

cr = {l,l, 1 , 1 , 1 ,- 1 ,- 1 , 1 , 1 ,- 1 , 1 ,- 1 , 1 }. 

The autocorrelation of x is given by 

A, = {1,0,1,0,1,0,1,0,1,0,1,0,13,0,1,0,1,0,1,0,1,0,1,0,1}. 

One can observe the peak and sidelobes of this autocorrelation in Figure 

The goal in radar coding theory is to construct code with zero sidelobes; however, in the case of Barker 
codes, this is impossible since the calculation of the last sidelobe, \Ax{N — 1 )| = |aiaA| = 1 , shows that A^ 
must always equal 1. On the other hand, if we employ a set of two complementary codes, i.e., a Golay pair, 
then it is possible to obtain zero sidelobes. This is illustrated in the next example. 
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Figure 1: Autocorrelation of Barker sequence x = {1,1,1,1,1, —1, —1,1,1, —1,1, —1,1}. 


Example 2.3. Consider two codes of length 4, xi = {1,1, — 1 , 1} and X 2 = {1, 1 , 1 , —!}■ Their autocorrela¬ 
tions are given by = { 1 , 0 , —1,4, —1, 0 , 1 } and = {—1, 0 , 1 , 4,1, 0 , — 1 }. Then {xi,X 2 ) forms a Golay 
pair since the sum of their autocorrelations has zero sidelobes: 


A,, + A,, = {0,0,0,4,0,0,01 


The notion of a Golary pair can be generalized to a set containing an arbitrary number of complementary 
codes (see |3]). 

Definition 2.2. A set of K codes {xi ,..., xk}, each of length N, is called a complementary code set if 

K 

Y,A.Aj) = NKS, ( 2 ) 

for —N -\-l<j<N—1 and where Sj is the Kronecker delta function defined by 


0, if,VO; 

' \l, ifj = 0. 

Next, we review Coxson and Haloupek’s [T] matrix formulation of complementary code sets. 

Definition 2.3. An NxK matrix M is a complementary code matrix (CCM) if its columns, Xi,..., Xjq, 
form a complementary code set. In that case, we define the composite autocorrelation function of M 
by the following 

K 

AM(j) = ^ A,,,(j) (3) 

SO that 

AmU) = NKS, (4) 

because of ©• 

Definition 2.4. Let M be a matrix. Then M ■ M*, where M* represents the conjugate transpose of M, is 
defined to be the row Gramian of Af. 
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Example 2.4. Let M = 
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-1 


. Then the row Gramian of M is given by 
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Definition 2.5. A N x N matrix Q is diagonally regular if its off-diagonals, i.e., those diagonals outside 
of the main diagonal, sum to zero. 


Note that Example 2.3 gives an example of a matrix M whose row Gramian is diagonally regular. 
Moreover, M is an example of a complementary code set. This not a coincidence. 

Lemma 2.1 (Goxson-Haloupek j^). An NxK matrix M is a CCM if its row Gramian is diagonally regular, 
i.e., 


M ■ M* = Q, 


where Q is diagonally regular. 

Next, we present the row-correlation function, which gives an equivalent representation of a GGM in 
terms of its rows [HIS]. 

Definition 2.6. The row-correlation function of a matrix M consisting of rows {ri,... jTk} is defined 
by 

N-j 

R-mO') = '^n- Ti+j 

i=l 

for 0 < j < N — 1 and RmU) = R-m(— j) for —N -I- 1 < j < 0. 

Note that RM^j) represents the sum of the elements in the j-th diagonal of M ■ M* , where j = 0 
corresponds to the main diagonal and j > 0 (or j < 0) corresponds to j-th diagonal above (or below) the 
main diagonal, respectively. Moreover, these row-correlation functions can also be represented as the sum of 
the autocorrelation function of each column in M. 

Example 2.5. The row-correlation function of M given in Example |2.4| is given by the following. 

Rm(0) =ri-fi+r2-f2 + rz-f3 

= [-1, -1, -*] • [-1, -1, *] + [-i, -i, 1] • [z, z, 1] + [-1, -1, 1] • [-1, -1,1] = [3,3, 3] 

R-m(I) = Cl • f2 -k r2 • fa = [-1, -1, -i] ■ [i, i, 1] -k [-i, -i, 1] • [-1, -1,1] = [-3i, 1 -f 2i] 

Rm(2) =ri-f3 = [-1, -1, -i] ■ [-1, -1, 1] = [2 - i] 

Rm(-1) =Rm( 1) = [3z,l-2z] 

Rm(-2) =Rm(2) = [2 + i] 

The following theorem characterizes GGMs in terms of the row-correlation function. 

Theorem 2.2 ([21 [B]). A N x K matrix M is a CCM if and only if 

K 

RM(j)=^A,.(j)=iVLfJ, 

i=l 

for —N + 1 < j < N — 1. 


4 







3 Symmetries of CCMs 

Let Cn,k{p) denote the set of all CCMs of dimension N x K consisting of p-th roots of unity. Any function 
/ : Cn,k{p) —t Cn,k{p) that maps CCMs to CCMs of the same dimension will be called a CCM-preserving 
symmetry of Ci<i,k(p)- In particular, if AmU) represents the composite autocorrelation of M, then we 
require 

^/(M)(j) = 0 

for all M e Cn,k{p) and |j| = 0,1 ,..., TV — 1. 

Coxson and Haloupek discussed the different ways to create a new N x K CCM using its symmetries as 
described in the following theorem [T]. 

Theorem 3.1. Suppose that M is ap-phase N xK CCM. Then the following operations are CCM-preserving 
symmetries, i.e., the resulting matrix is a N x K CCM. 

(i) Take any eolumn x of M and replaee, x, with ax where a is any unit-modulus complex number. 

(ii) Take any eolumn x = [mi,..., m^] of M and replace x with its conjugate reversal x = [rriN-, ■ ■ ■, fhi]. 

(Hi) Form M. 

(iv) Form the product Diag\l3, fP, 0^,... with /3 being any unit-modulus complex number. 

(v) Form MP where P any K x K permutation matrix. 


Using this theorem we can see that with one CCM many more can be created. Coxson and Russo created 
an algorithm for narrowing down the search space for complementary code matrices using these symmetries 
[5] . This is important due to the fact that the search space for p-phase N x K CCMs is p^^] as the 
dimensions increase the search takes exponentially longer to conduct an exhaustive search. For example, 
the exhaustive search for 4-phase 4x4 CCMs, which generated 4300800, took about 24 hours using parallel 
processing on 4 cpus. With the help of this theorem, it is possible to narrow down the search space if it is 
possible to advoid checking all but one from a set of equivalent CCMs. There is a wonderful example of this 
by Gibson where he classified 4-phase Golay pairs into equivalence classes [5]. We aim to do the same, but 
for more general CCMs. 

A quick example of this is as follows. 


Example 3.1. You can see here that we can take one matrix M and make two other matrices out of M, 
though we can make many more besides these two using Theorem 3.1 To get the first matrix we applied (i) 
multiplying xi,X 2 , and X 4 by —i and X 3 by —1. Then to that matrix we applied (iv) with /3 = —i. Yet the 
other interesting thing about this example is that we were able to make the entire first row ones as well as 
777 . 2 , 1 . This is what we call a normalized CCM. 
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3.1 Complementary Codes 

We proceeded to find the relationships between different CCM symmetries. Our results generalize Coxson’s 
who examined symmetries that preserve low sidelobe levels in Barker codes [7]. We begin with definitions 
and notations to express the symmetries mentioned in Theorem |3.1| 

Definition 3.1. Let the following assist in representing the different operations Theorem 3.1 

(a) Z = {( 01 , 02 ,...,Oif) : = 1 }. 

(b) A = ^3, ■ • ■, tx) ■■ tk e {0,1 }}. 

Definition 3.2. The following are definitions that represent the CCM preserving operations. Theorem |3.1| 
where U G Z and T G A. 

(a) C denotes multiplying columns in M Theorem |3.1| (i) 


CuM = [oicci, 02X2,..., ukXk] 


(b) p denote the conjugate reversals of columns in M. Theorem 3.1 (ii) 


PtM=[u}‘^ {X 2 ), ■ ■ ■, 


(c) S denote the conjugate of the matrix M. Theorem 3.1 (iv) 


SM = M 


(d) Q{P) = N X N diagonal matrix with values going from /3, /3^,... Theorem |3.l| (iv) 
Q{I3)M = [Pri, 13^X2,...,= [Q{l3)xi,Q{l3)x2, ■ ■ ■ ,Q{P)xk] 


(e) P = a permutation of the columns. Theorem 3.1 (v) 


P^ ; Xfj^2) ; ^(t(3) 7 • ■ • ; 


(f) t/ = [ai, 02 , 03 , . . . , Olc] 

(§) Up 7 ^(t( 2 ) 7 ^ct( 3 ) 7 ■ ■ ■ 7 ^cr(/C)] 

(h) Ut,p = [oi, 02 7 Q!3, ..., Oix] such that where tk G T 




We now establish relationships between symmetries by viewing them as group generators. 

Lemma 3.2. Let M be a N x K p-phase CCM, M = [xi,X 2 , X 3 ,..., xk] and 
Xk = [mi^k,m 2 ,k,'m 3 ^k, ■ ■ ■,mN,k]'^. Then 

(i) CjjPtM = prCurpM 

(li) CuSM = SCjjM 
(ill) CuQW)M = Q{P)CuM 
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(iv) CuPM = PCupM 

(v) ptSM = SpxM 

(vi) ptPM = PpTp_^M 

(vii) SQ{P)M = Q{P)SM 

(viii) SPM = PSM 


(ix) Q{(3 )ptM = Cut^uPtQpM 

(x) Q{fi)PM = PQ{I3)M 
proof. (i) 

CuPtM =CuPt[xi,X2,X3,...,xk] 

= Cu (xi), (x^), a;2-*3 (^ 3 ),..., ( 3 ;^)] 

= pT[uj'^~*'Hai)xi,uj‘^~'^'^ {a2)x2,uj‘^~^^ {a3)x3,... {aK)xK\ 

= PtC'c/^ [xi,X2,X3, . . . ,Xk] 

= PtCjJt^ 


(v) 


PtSM = pt-S'[xi,X 2 ,X 3 , ... jX/f] 

= Pt[xi,X2,Xz, ■ . . ,Xif] 

= (xi), (^ 2 ), C.2-‘3 (^ 3 ), . . . , (^^)] 

= (xi), (x2), w 2-‘3 (3,3)^ . . . ^ ^2-t^ (^^)] 

= Spt[Xi,X2,X3, . . .,Xk] 

— SptM 


(vi) 


prpPM = PtP[xi,X2,X3,...,Xk\ 

PT\Xa{l) j Xfjp2^ , a7(T{3) ? • • ■ 5 

= [w2-*i(x„(i)), w2-*3(x„(2)), w2-*3(x^( 3)), ..., (a:<T(K))] 

= w^~*‘'“^( 3)(2;3)^ _ ^ (x*:)] 

= PpTp.i [xi,X2,X3, ...,Xk\ 

= PpTp.^M 


(ix) 


Q(/3)ptM = (3(/3)pT[a:i,X2,X3,... ,xk] 

= Q(/3) (xi), w 2-‘2 ( 3 ^ 2 ), a; 2-*3 ( 0 : 3 ),..., {xk)] 

= [Q(^)w^"‘' (a;i), Q{P)uj‘^~*^ (X 2 ), (3(^)w2-*3(X 3 ),..., Q{l3)uj‘^-*^( xk)] 

= Cut^p [aiQ(^)w^“‘i (xi), q; 2 Q(/ 3 )w^“‘= (X 2 ), a 3 ( 5 (/ 3 )a; 2 “* 3 ( 3 - 3 )^ _ ^ aKQ{P)oo'^~*‘^ (xk)] 

= C';7^_^pr[w^“‘3(Q;i(5(/3)a;2“*3(xi)),a;2“*3(a2(5(/3)w^“*3('3,2))^, _., (a_R-Q(^)w 2 “‘^(x_r-))] 

= C';7^^pr[w^“‘3(Q;i(5(/3))xi,a;2“*3('^2(5(^))x2,.. . {aKQ{f3))xK] 

= Cu^ i,PtQ{P)[xi,X2,X3, . . . ,Xk] 

= Cut^isPtQpM 


The proofs for parts (ii), (hi), (iv), (vii), (viii), and (x) are analogous and will be omitted. □ 

Definition 3.3. The complementary group G of the set of all x it! p-phase CCMs is defined to be the 


group generated by the symmetries S, P,Cu, PT,Qi/3) and their relations given in Lemma 3.2 
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Theorem 3.3. The cardinality of the complementary group G of Cn,k{p) is bounded by 


\G\ < 2^+^p^+^K\ 


(8) 


Proof. Each CCM preserving operation on the matrix M can be represented in the following form. 

SPGuPtQW) 

Since S represents the conjugate of the matrix there are two possibilities, i.e., IS”! = 2. P represents 
permutations of the columns |P| = K\. Cu is multiplying columns by different unit-modulus complex 
number i.e. \Gu\ = P^■ Whether or not the conjugate reversal of a column is take has \pt\ = 2,^■ And if we 
do progressive multiplication there is |Q(/3)| = p- So this will will produce a max of 

|G| < \S\\P\\Gu\\pT\\Qm = 2K\p^2^p = 2^+^p^+^K\ 


CCMs from M. 


□ 


4 CCM Search Algorithm 


The following is an explanation of our exhaustive search algorithm for CCMs. Our search algorithm is based 
on the search algorithm of Coxson and Russo |2]. While their work searches for the binary CCMs, ours 
focuses on any poly-phase. We begin by presenting how a CCM can be represented in a normalized form as 
shown in the following lemma an example of which was expressed in Example 


Lemma 4.1 (Normalized CCM). Any N x K CCM can be normalized by transforming the matrix into 
another matrix (using Theorem 3.1) where the entire first row consists of Is and the first element in the 
second row is a 1. 


Proof. Let M he a, N x K CCM with ri = [mi i, mi 2 , mi_ 3 ,..., mi^x] and r 2 = [m 2 ,i, m 2 , 2 : m 2 , 3 ,..., m 2 ,if]. 
In this proof it is sufficient to just show the operations being applied to the first two rows, but the reader 
may note that these operations will be applied to the other rows as well. 

Here we used (i) in Theorem |3^ multiplying each column of M by the conjugate of the first value of each 
column. This makes the entire first row ones since xx = 1. 


( mi,i mi ,2 mi ,3 

m2,i m 2,2 m2,3 


mi^K 

m2,K 


/ mi,imi,i 
mi,im2,i 


mi,2mi,2 mi,3mi,3 
mi,2m2,2 mi,3m2,3 


V 


mi,Km2,K 


Next we uses (iv) in Theorem 3.1 and the fact that Diag(m 2 ,imi,i) =Diag(m 2 ,imi,i). 


( m2,imi,i 

(m2,imi,i)^mi,im2,i 


m2,imi,i m2,imi,i 

(m2,imi,i)2mi,2m2,2 (m2,imi,i)2mi,3m2,3 


m2,imi,i 

(m2,imi,i)2mi,Km2,if 


The final step is to now multiply each column by m 2 ,imi,i = m 2 ,imi,i to bring the first row back to all I’s. 
Note now though that in this new matrix, M', m^ 1 = (m 2 ,imi,i)^(mi,im 2 ,i)^ also known as m 2 1 = 1 The 
normalized version of M can now be written as M'. 


M' = 


( 1 


V : 


1 

’7^2,2 


^2,3 


m 


2,K 


□ 








4.1 Short summary of search algorithm 


The algorithm is an exhaustive search for all representations of CCMs. The search implements the known 
symmetries to avoid finding a mass amount of equivalent CCMs. We construct the matrix from the outside 
in. We start with making sure the matrix is normalized, setting the first row all equal to ones and first 
element of second row a one. We then force the last row to sum to zero. Next we imagine the unit circle and 
name the first angle, between the origin and the first of the p-phase, uj. We can define a method of sorting 
where Out < uj < 2uj <...< {p — l)a;. When constructing r 2 we want the corresponding if-tuples to follow 
this sorting process. We then construct r^v-i by only using if-tuples whose sum is equal to the opposite of 
the partial sum of the row-correlation function of the second to last off diagonal, meaning summing up all 
of the off diagonal minus the one that includes the row that we are on. We save time here by not looking 
at if-tuples that we know will make the — 2 off diagonal not sum to zero. From here we create based 
off of r 2 . If r 2 has consecutive values that are equal, then at those positions in the values need to be 
sorted. Finally we continue searching for values for the other rows by continuously looking at their partial 
sum of their row-correlation function. Once we have created all of these rows we just need to check that the 
remaining TV — 1 — y off diagonals also sum to zero. If so the matrix is a COM. 

There are many similarities and one major difference between our search algorithm and Jon Russo’s. The 
similarities consist of constructing the code from the outside in and calculating the autocorrelation functions 


as the search progresses. Both algorithms normalize the first row, to account for Theorem 3.1 (i), and set 
the sum of the last row to zero. They also make sure the conjugate reversal of the column is less than 


the column, accounting for Theorem 3.1 (ii). Both also make sure that the columns are ordered to avoid 


permuting the rows. The difference between the search algorithms is that ours is generalized for all N x K 
poly-phase CCMs. 

One of the most difficult parts of this algorithm to understand is the sorting by increasing angles. For 
this reason we will begin with an example for the reader. Another way one might think of the increasing 
angles is by increasing exponents. In a 4-phase CCM we deal with {1, i, —1, —i} or {i°, z^}. 

Definition 4.1. Let r„ be the n-th row of a matrix M. 

Example 4.1. Let r„ = [1,1, z] Then the options for the increasing exponents of r„ are as follows: 
{[1,1,1], [l,z,l], [1,-1,1], [l,-i,l], [i,f,l], [z,-l,l], [f,-f,l], [-1,-1,1], [-!,-*,!], [l,lo], 

[l,b*], [1,-1,*], [1,-«,*], [l*a], [b-1,*], [-l,-lo], [-1,-b*], [-b-b*], [1,1,-1], [l,b-l], 

[1,-1,-1], [l,-z,-l], [z,z,-l], [z,-l,-l], [z,-z,-l], [-1,-1,-1], [-l,-z,-l], [-z,-z,-l], [1,1,-z], 

[l,b-i], [1,-1,-*], [1,-*,-*] [*,*,-*], [*,-!,-*], [*,-*,-*], [-1,-1,-*], [-1,-*,-*], [-*,-*,-*]} 

Instead of there being 4^ possibilities for the next row there are 3x4^. As the rows get longer or r„ 


changes, this can save more and more time. All we are doing here is sorting the columns as Theorem 3.1 (v) 
allows. 

The beauty of constructing the matrix from the outside in is that this construction method allows for 
the next key part of this algorithm. We can sum up the elements off diagonals of the partially constructed 
matrices dot product with its conjugate transpose. If we sum up all of the terms except the term that is 
multiplied to row one. This sum will first have to be less than or equal to K by in order for this to happen. 


( r^\ 
r2 
rs 
V *-4 y 


ri r2 rs ) = 


/ rifi 

rif2 

*■1^3 

rif4 \ 

r2ri 

*•2^2 

*■2^3 

*■2^4 


*’3^2 

*■3^3 

rsTi 

\ *’4*’l 

*’4*'2 

*■4*^3 

r4f4 / 


(a) Increasing exponents (E(r„)): creates a list of new rows such that for each row r ^', = [z*'^, ■ ■ •, *^^] 

such that is sorted based off of r„. Represent r„ = [to„_i, to„^ 2 , ■ ■ •, **Zn,if] for any consecutive equal 
m’s in the r„' of those values must be in increasing exponential form mod 4. If happens to be 
normalized then E(r„) = a list of r^Fs where = [z^’L*^^, • ■ • ,*^^] and mod 4 < Pk+i mod 4. If 
r„ = [1, l,z, —z] = [z°,z°,z^,z^] then E(r„) = a list of r^/’s where r^' = [z*’L*^’^, • ■ • ,*^^] where pi mod 
4 < P 2 mod 4 and 0 < pa < 3, 0 < p 4 < 3. 
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(b) Define the absolute value of a complex number as the Taxicab Distance: ||(a + 6i)|| = |a| + |6| 


4.2 The Search 

This search was implemented in C++ using the adaptation of the Coxson-Russo search algorithm mentioned 
in this section. Since we used the symmetries to speed up the search, we did not exhautively apply the 
symmetries. So once the search was completed, a Mathematica program was implemented on the matrices 
found where all combinations of symmetries being applied to them. The purpose was to narrow it down to 
the exact number of equivalence classes. The results of which are shown in Table [l] 

5 Equivalence Classes of CCMs 

Definition 5.1 (Equivalent CCMs). Let Ml, M 2 S Cn,k{p) and g G G, where G is the complementary 
group. Then Mi and M 2 are said to be equivalent if 

g{Mi) = M 2 

We define Cm to be the equivalence class of Gm,k{p) containing M, i.e.. 

Cm = {M G Cn,k{p) ■ gM = M,g G G} 

The number of V x 2 4-phase CCMs for N = 1, 2, 3,..., 22 are already known as previously discussed [3]. 
So we will look at what this definition means for N x A CCMs. The order of the symmetry group of 4-phase 
V X 4 CCMs 786432. We must check each of the resulting matrices to ensure that the set of CCMs is 
narrowed down to the simplest number of equivalence classes. With this theorem we able to narrow down 
the 4-phase 4x4 CCMs to 24 and 4-phase 6x4 CCMs to 1448 with the assurance that these are the number 
of equivalence classes, no more, no less. We also found the number of equivalence classes for other, easier 
N’s. We then looked at the equivalence classes to see how many of them could be formed using a Hadamard 
matrix. It is interesting to mention here that the entire equivalence class need not be Hadamard matrix in 
order to have a Hadamard matrix representation. This is because taking the conjugate reversal of a column 
does not always preserve the Hadmard matrix property 


Table 1: Equivalence Classes 


N X K p-phase CCMs 

*CR Algorithm 

Number of Equivalency Classes 

Hadamard Representations 

2x4 4-phase 

36 

2 

2 

3x4 4-phase 

95 

5 

5 

4x4 4-phase 

231 

24 

17 

5x4 4-phase 

5246 

133 

0 

6x4 4-phase 

23448 

1448 

0 


*Coxson and Russo’s adapted search algorithm 

5.1 Constructing Quaternary CCMs 

When dealing with CCMs, many times it is useful to look at constructions that can be used to create larger 
CCMs from smaller ones since the search space for the smaller CCM is, of course, smaller. It is here that the 
reader could think back to some of the more well known constructions. The Kronecker product, represented 
by the symbol 0, can take matrix Mi that is a W x Ki CCM and M 2 that is a iV 2 x K 2 CCM to create 
a new matrix M where Mi 0 M 2 = N 1 N 2 x KiK 2 . In the concatenation theorem you can take Mi to 
he & N X Ki CCM and M 2 to be a iV x K 2 CCM. The concatenation theorem gives us the following: 
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Algorithm 1 Search Algorithm for N x K CCMs 

Normalize ri {from this point forward we will write fi as ri} 
for all K-tuples that sum to zero do 
set tn equal to this K-tuple 

if A = 2 then 
This is a CCM 
end if 

Set 7712,1 = 1 {now the matrix is fully normalized} 
for all K-tuples in E(ri) do 
set r 2 equal to this K-tuple; {Partial sum of the rows 2} 

if A = 3 then 

if i?M(l) = 0 then 
This is a CCM 
else 

NOT A CCM 

end if 
end if 

if ||Rm(2 - A) - rjjv-qT-ill < K then 

for all K-tuples that are equivalent to — (Rm( 2 — A) — r[iv_i]ri) do 
set tn-i equal to this K-tuple 

while 2 < k < K, m 2 ,k = and mN-i,k''nM,k = where 0 < P < 3 and 0 < P' < 3 do 

if P > P' then 

Not the CCM we are looking for; {This step is taking into account the conjugate reversals 
of a column} 

end if 
end while 

for all K-tuples in E(r 2 ) do 
set ra equal to this K-tuple 
for t = l,t < — 2, t + P do 

if ||RM(f + 2- N)- r[Ar_t_i]ri|| < K then 

for all K-tuples whose sums are equivalent to —(RM(t + 2 — A) — r[Ar_t_i]ri) do 
set equal to this K-tuple 

end for 
end if 
end for 

for s = l, s<A—1—^,s-|--|- do 
if Rm(s) 7 ^ 0 then 

NOT A CCM; {Check that the remaining off diagonals sum to zero} 

end if 
end for 
This is a CCM 
end for 
end for 
end if 
end for 
end for 
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= Nx {K 1 +K 2 ) H]. There are even more that you can read about yet all the ones from Coxson 
and Haloupek focus, as we said, on creating a CCM from a smaller dimension CCM to decrease the search 
space [1]. We can also look at creating a quad-phase CCM from the ternary CCMs. These are a looser 
definition of CCMs whose entries consist of the values {—1, 0,1} and their off diagonals sum to zero. Here 
we decrease the search space, 3^^ as opposed to 4^^. Say we have a matrix M that is a 4-phase CCM. 

/ -1 -1 -i -i \ 

M = * * 

i -i -1 1 

\ i -i i -i J 

We can decompose M into the matrix’s real and imaginary parts M = A + iB, where A and B are ternary 
CCM’s. 



-1 

0 

0 \ 


f ^ 

0 

-1 

\ 

0 

0 

0 

0 

,B = 

-1 

-1 

1 

1 

0 

0 

-1 

1 

1 

-1 

0 

0 

V 0 

0 

0 

0 ^ 


1 1 

-1 

1 

-1 / 


Though using the ternary CCMs may seem trivial, they can have power. 

Theorem 5.1 (Dual Pair Theorem). Assume that A and B are NK ternary CCMs. Then Z = A + iB 
is a N X K quad-phase CCM if 

(i) I An.k I + I Bn k 1= 1 {Vn, k\l <n < N and 1 <k < K} 

(li) BA* — B*A is Diagonally regular 

NOTE: Any pair of matrices (A, B) that satisfy condition (i) is called a dual pair. 

Proof. We first prove that Z is quad-phase by considering two cases: 

Case I'.Aij = 0 and Bi j = ±1. Then Zij = Aij -\- iBij ^ Zij = 0 ± i ^ Zij G {1, i, —1, —i}. Hence, Z is 
quad-phase. 

Case 2: Ai j = ±1 and Bi j = 0. Then Zi j = j -|- iBi j ^ Zi j = ±1 -|- iO ^ Zi j G {1, i, —1, —*}. Again, 
this proves that Z is quad-phase. 

Next, we prove that ZZ* is diagonally regular by first calculating 

ZZ* = {A + iB){A + iB)* 

= {A + iB){A* -iB*) 

= AA* -h i{BA* - B*A) + BB* 

Since A and B are ternary CCMs, AA* and BB* are diagonally regular. Moreover, by assumption BA* — B*A 
is also diagonally regular. Thus, ZZ* must be diagonally regular. □ 

In the following table, a large portion, and in some cases all, of our CCM equivalence classes can be 
represented as a dual pair, which has a smaller search space. This is represented in the following table that 
compares the results of implementing different types of construction methods. 


Table 2: CCM Equivalence Classes 


N X K p-phase CCMs 

Equivalence Classes 

A© H 

[A,H] 

Dual Pair Representatives 

2x4 4-phase 

2 

n/a 

2 

2 

3x4 4-phase 

5 

n/a 

1 

5 

4x4 4-phase 

24 

2 

6 

22 

5x4 4-phase 

133 

n/a 

3 

94 

6x4 4-phase 

1448 

2 

27 

471 
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7 Appendix 

I. 2 X 4 Equivalence Class Representations: 

1 . [[ 1 , 1 , 1 , 1 ], [ 1 , 1 ,- 1 ,- 1 ]] 

2. [[l,l,l,l],[l,i,-l,-i]] 

NOTE: Both classes can be represented by dual pairs. 

II. 3x4 Equivalence Class Representations: 

1 . [[ 1 , 1 , 1 , 1 ], [ 1 ,- 1 ,- 1 ,- 1 ], [ 1 , 1 ,- 1 ,- 1 ]] 

2. [[l,l,l,l],[l,-l,-l,-l],[l,i,-l,-i]] 

3. [[l,l,l,l],[l,i,-l,-i], [1,-1,1,-1]] 

4. [[l,l,l,l],[l,-l,-l,-i],[i,i,-i,-i]] 

5. [[l,l,l,l],[l,i,i,-!],[-!, 1,1,-!]] 

NOTE: All classes can be represented by dual pairs. 

III. 4x4 Equivalence Class Representations: 


1. 

[[1.1,1,1], [1,1,1 

,1],[1,1,-1,-1],[-1,- 

-1,1,1]] 

13. 

[[1,1,1,1], [1,1. 

*,-1],[1,-*, - 

1,1].[*.- 

-1,1.-i]] 

2. 

[[1.1,1,1], [1,1,1 

,i],[-l,-l,-i, !],[!,- 

1,1,-!]] 

14. 

[[1,1,1,1], [1,1. 

*,-1],[1,-*, - 

*,*],[*,- 

-1.1.-i]] 

3. 

[[1.1,1,1], [1,1,*, 

*],[i,-i,b-*],[-i,i. 

1,-1]] 

15. 

[[1,1,1,1], [1,1. 

1,*], [*,-!,-*, 

-*],[-*: 

,bi.-i]] 

4. 

[[1.1,1,1], [1,1,*, 

*],[*,-*, 1,-1], [1,-1, 

-1,1]] 

16. 

[[1,1,1,1], [1,1. 

l.*],[-l.-l.- 

-*. !],[*,- 

-bi.-i]] 

5. 

[[1.1,1,1], [1,1,- 

-i,-i],[i,-i,i,-i],[i 

,-l,-l,l]] 

17. 

[[1,1,1,1], [1,1. 

*, -i],[b-i,i, 

-1],[1, 


6. 

[[1.1,1,1],[1,*,*, 

-!],[-!,i,i, 1],[-1,1, 

1,-1]] 

18. 

[[1,1,1,1], [1,1. 

*.*], [*,-*,*.-> 

;],[-*,*, 

l.-l]] 

7. 

[[1,1,1,1], [1,1,1 

,l].[l,i,-l,-i],[-l,- 

i, l,i]] 

19. 

[[1,1,1,1], [1,1. 

1,1], [*,*,-*,- 

*], [-*, - 

-*,*,*]] 

8. 

[[1,1,1,1], [1,1.*, 

i], [1,-1,1,-1], [-1,1, 

-b*]] 

20. 

[[1,1,1,1], [1,1. 

-1,-1], [*,-*, 

*, -*],[*: 

,-*,-*,*]] 

9. 

[[1,1,1,1], [1,1,*, 

*],[*,-1,1,-*], [1,-1, 

-i,i]\ 

21. 

[[1,1,1,1],[1,*. 

*,-i],[-*,-i. 

-!.*],[- 

-*,*,*,-*]] 

10. 

[[1,1,1,1], [1,1.*, 

*],[-!,-1,1,1], [1,-1, 

b-*]] 

22. 

[[1,1,1,1],[1,*. 

-!,-*],[-*,*,- 

-*,*], [- 

*,-l,*,l]] 

11. 

[[1,1,1,1], [1,1,*, 

*],[*, -b*, -*],[i,-i,- 

-i,i]] 

23. 

[[1,1,1,1],[1,*. 

-l.-*],[-l,l. 

-1,1], [ 

-l,i, 1, —* 

12. 

[[1,1,1,1], [1,1,- 

-i,-i],[i,-i,i, -i],[i, 

-1,-*,*]] 

24. 

[[1,1,1,1],[1,*. 

-l.-*],[l,-l. 



NOTE: All classes, except classes 23 and 24, can be represented by dual pairs. 

IV. 5x4 and 6x4 Equivalence Class Representations: These representations can be downloaded from our 
website at the following address: http://elvis.rowan.edu/datamining/ccm/equivalence/ 
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